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ABSTRACT 

The next generation of long baseline neutrino oscillation experiments will aim at 
determining the unknown mixing angle ^13 , the type of neutrino mass hierarchy 
and CP- violation. We discuss the separation of these properties by means of the 
energy dependence of the oscillation probability and we consider an hybrid setup 
which combines the electron capture and the decay from the same radioactive 
ion with the same boost. We study the sensitivity to the mixing angle and the 
CP-phase, the CP discovery potential and the reach to determine the type of 
neutrino mass hierarchy. The analysis is performed for different boosts and 
baselines. We conclude that the combination of the two decay channels, with 
different neutrino energies, achieves remarkable results. 

1. What is known, what is unknown 

In the past d ecade, atmospheric , solar l^MSJ ^ reactor l ^ l '^ l ^ l' and long-baseline 
accelerator I^SUSI neutrino experiments have provided compelling evidence for the phe- 
nomenon of neutrino oscillations. This has reshaped our understanding of the prop- 
erties of elementary particles as it implies that neutrinos have mass and mix. The 
combined data can be described by two mass squared differences, Amg^^ and Am^i, 
where Am^j = — m^, whose current best fit values are [Amg^^l = 2.4 x 10"'^ eV^ 

and Am2i = 7.65 x 10~^ eV^^^^. The two mixing angles 612 and 621, drive the solar 
and KamLAND, and atmospheric and accelerator neutrino oscillations, respectively, 
and are measured to be sin^ ^12 = 0.304 and sin^ ^23 = 0.50^^. The third mixing 
angle, 6'i3, is yet undetermined but is known to be small or zero. With available data, 
6^13 is constrained to be 

sin^ ^13 < 0.040 (0.056) at 2(t (3(t) . (1) 

It is interesting to note that very recently a first hint in favour of ^13 7^ has been 
found in a combined analysis of atmospheric, solar and long-baseline reactor 
neutrino data, with: 

sin^ ^13 = 0.016 ±0.010 at la , (2) 

implying a preference for 6*13 > at 90% CL. The combination of ([2]) with other 
analyses is discussed by G. L. Fogli in the present Proceedings. 



Although the experimental progress in neutrino physics over the last decade has 
been conspicuous, many of the fundamental questions surrounding neutrinos still need 
to be addressed. Understanding of the physics beyond the Standard Model responsible 
for neutrino masses and mixing requires knowledge of the nature of neutrinos (whether 
Dirac or Majorana particles), the neutrino mass ordering (normal or inverted), the 
absolute neutrino mass scale, the value of the unknown mixing angle ^13, and whether 
CP-symmetry is violated in the lepton sector. It will also be necessary to improve 
the precision on the known parameters, in particular to measure any deviation from 
maximal ^23 and, if so, to determine its octant. 

Some of the iss ues ab ove will be addressed by a future program of neutrino oscil- 
lation experiments l^^^l^. For flavour oscillations, the unitary diagonalization of the 
neutrino mass matrix, assuming that the flavour mixing affects the three active light 
neutrinos only, is given by the PMNS matrix U which connects from mass eigenstate 
neutrinos to flavour eigenstate neutrinos 
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U is determined by 3 mixing angles and 1 CP-phase, even if neutrinos are Majorana 
particles. The two additional Majorana phases would need a A(L) = 2 process to 
become observable, as in Neutrinoless Double Beta Decay. 

Nuclear reactors EHl ^nd long baseline experiments using conventional beams™ 
will be the first to explore ^13 below the current limit and maybe confirm the hint for 
6^13 7^ . If 6*13 is close to the present bound ir nposed by the running and near future 
experiments, the next generation of superbeams , an extension of a conventional 
beam with an upgrade in intensity and detector size, and wide-band beams will 
probe CP- violation and, for sufficiently long baseline, the neutrino mass hierarchy. For 
small values of 613 or, if ^13 is large but a better precision on the neutrino parameters 
needs t o be a chieved, the connnunity must turn to the novel concepts of the neutrino 
factory l 30 | 21 | qj. beta-beam I 22 | 23 | whereas conventional beams sourced from pion 
decays have an intrinsic contamination of electron neutrino at the ~ 1% level (owing 
to kaons in the beam), neutrino factories and beta-beams will have clean sources from 
highly accelerated muons and ions, respectively, producing a well-collimated beam. 
In a neutrino factory, muons (antimuons) are produced, cooled and accelerated to a 
high boost before being stored in a decay ring. The subsequent decay sources a muon 
neutrino (muon antineutrino) and electron antineutrino (electron neutrino) which are 
aimed at magnetised detectors located a very long distance from the source. The use 



of magnetised detectors is necessary to separate the 'right muon' disappearance signal 
from the 'wrong muon' appearance signal, which is sensitive to matter effects and CP- 
violation. A beta-beam will exploit accelerated ions that (3-decaj sourcing a clean, 
collimated, electron neutrino beam. Magnetised detectors will not be necessary in 
this case, the only requirement being possession of good muon identification to detect 
the appearance channels. Therefore, water Cerenkov (WC), totally active scintillator, 
liquid argon detectors and non-magnetised iron calorimeters could be used, depending 
on the peak energy. 

The determination of the oscillation parameters is severely affected by degenera- 
cies: the possibility that different sets of the unknown parameters {sgn{Amli) , 6, 6*13, 623 
octant) can provide an equally good fit to the probability for neutrino and antineutrino 
oscillations, for fixed baselines and energy. Therefore, a high precision measurement 
of the appearance probabilities is not sufficient to discriminate the various allowed 
solutions. In order to weaken or resolve this issue, various strategies have been put 
forward: exploiting the energy dependence of the signal in the same experiment, us- 
ing reactor neutrino experiments with an intermediate baseline, combining different 
long baseline experiments, adding the information on ^13 from reactor experiments, or 
using more than one baseline for the same beam. In addition, ^13 controls the Earth 
matter effects in multi-GeV atmospheric and accelerator neutrino oscillations, as well 
as in supernovae. These might provide useful information on the type of neutrino 
mass hierarchy and ^13. 

2. CP- violation in neutrino oscillations 

The magnitude of the T-violating and CP-violatin g inter ference in neutrino os- 
cillation probabilities is directly proportional to sin ^^3^2412511 CP- violation can be 
observed either by an Asymmetry between neutrinos and antineutrinos and/or by 
Energy Dependence in the neutrino channel. In the last case, the CP phase 6 plays 
the role of a phase shift in the interference pattern between the atmospheric and 
solar amplitudes for the appearance oscillation probability. This result is a conse- 
quence ^^Si Qf ^]-^g assumptions of CPT-invariance and No Absorptive part in the 
oscillation amplitude: the Hermitian character of the Hamiltonian responsible of the 
time evolution says that the CP-odd=T-odd probability P(z/e — > i^fj) — -P(t'e ^v) is 
an odd function of time, i.e., an odd function of the baseline L. In vacuum neutrino 
oscillations for relativistic neutrinos, the oscillation phase depends on the ratio L/E, 
then the CP-odd term becomes an odd function of the energy E for fixed L. With 
the same reasoning, the CP-even terms are even functions of the energy E in the 
oscillation probability. In this way. Energy Dependence in the appearance oscillation 
probability is able to disentangle CP-even and CP-odd terms. 

The explicit expression for the suppressed appearance probability for neutrinos in 



vacuum oscillations is given by 
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where A12 = Aml^/{2E), A13 = Aml^/{2E) and J = C13 sin 26*12 sin 26123 sin 26*13 . 
This expression shows, in fact, that |[/(e3)| gives the strength of the probability, 
whereas S governs the interference pattern as a phase shift. Furthermore, the sep- 
arate atmospheric probability, the solar probability and the CP-even term of the 
interference are even functions of E/L. On the contrary, the CP-odd term of the in- 
terference between the atmospheric and solar amplitudes is an odd function of E/L. 
This result suggests the idea of disentangling 6 from |f/(e3)| without a need of com- 
paring neutrino and antineutrino events, which have different beam systematics and 
different cross sections in the detector: either monochromatic neutrino beams with 
different boosts or a combination of channels with different neutrino energies in the 
same boost are able of separating the CP-violating phase. 

Due to neutrino propagation through the Earth, matter effects can "fake" CP- 
violation in the sense that the presence of matter affects neutrino and antineutrino 
oscillations in a different way. It is not easy to disentangle matter effects from CP- 
violation since there is the so-called "mass hierarchy degeneracy", which swaps the 
effect of matter for neutrino and antineutrino oscillation according to the sign of 
Am|i. The energy dependence of the z/g Ufj, os cillatio n probability, in presence of 
matter effects, can be studied from the expression l ^^ ^ ^'^l ' 

Pi^e ~^ ^fi, L) = Tatm + Tsol + Tint , 



where A = ^J^GpneiJ-) and fiJ^V) = 1/L J^^ ne{L')dL' is the average electron number 
density. As seen, the energy dependence induced by the presence of A 7^ is different 
in Tatm, Tsol and in the interference Tint and, in fact, different from the energy depen- 
dence associated with the CP-even versus the CP-odd separation. On the other hand, 
the mass hierarchy degeneracy in vacuum is now removed because Tatm and T^ are 
now changing under the change of sign of Am|i, although Tsoi remains the same. All 
in all, we observe the virtues of studying the neutrino appearance probability as a 
function of the neutrino energy. 

From the current discovery phase of 6*13, next generation experiments will hence 
aim at precision measurements of the u^. — > z/^j oscillation probability. This will require 
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large underground detectors coupled to more intense and pure neutrino beams. These 
aspects are being studied within the LAGUNA and EURONu design studies. The 
knowledge on the possible values of ^13 is a necessary input to best optimize the search 
for CP-violation in the leptonic sector. 

The energy dependence of the signal is typically used to extract information on 
the mass hierarchy and CP-violation. Matter effects increase with baseline and en- 
ergy suggesting that setups with baselines L> 600 km are necessary I 28 | 29 | 30| ^^le 
determination of the type of neutrino mass ordering. In beta-beam experiments, such 
strategies would make use of a proposed upgrade of the CERN Super Proton Syn- 
chrotron (SPS) which would equip the accelerator with fast superconducting magnets 
allowing high boosts and fast ramps. The latter are important to reduce the loss of 
ions through decay in the acceleration stage. A sister approach to the beta-beam is 
to us e the neutri nos sourced from ions that decay mainly through electron capture 
(EC) I 31 | 32 | 33 [34li ^.j^g electron capture decay is dominated by a single channel, 
then a monoenergetic electron neutrino beam can be produced this way. In this case, 
all the beam intensity can be concentrated at the appropriate energy to get the best 
sensitivity to the oscillation parameters. In order to disentangle the CP violating 
phase with neutrino beams only, one makes use of the different energy dependence of 
the CP-even and CP-odd terms in the appearance probability Electron capture 
competes with P~^-decay when the QEC-value > 2me, rrie being the electron mass. 
With the ions identified in'^, the use of an upgraded SPS or the Tevatron0 is nec- 
essary to source baselines in excess of CERN-Frejus (130 km). We will discuss here 
also a hybrid of these two approaches. By selecting a nuclide with Qec ~ 4 MeV, 
we can make use of neutrinos from an electron capture spike and f3~^ continuous spec- 
trum simultaneously. Assuming a detector with low energy threshold, the use of such 
ions allows one to exploit the information from the first and second oscillation max- 
ima with a single beam. The use of the hybrid approach makes it possible to use 
a monochromatic beam at higher energies and a beta-beam at lower energies. The 
need for good neutrino energy resolution at the higher energies will therefore be less 
crucial than for high-7 beta-beam scenarios. 



3. A combined beta-beam and electron capture neutrino experiment 

In this section we discuss the 

idealSSl 

of the beta-beam and electron capture hybrid 
approach. We present the spectra of the two branches, their ratio and consider two 
nuclides which have desirable properties. 

The beta-beam is a proposal, originally put forward by P. Zucchelli to ac- 
celerate and then store /3-emitting ions, which subsequently decay to produce a well 
coUimated, uncontaminated, electron neutrino (or antineutrino) beam. The high lu- 



'^Note that the present Tevatron configuration does not ramp fast enough, resulting in a high loss 
of ions, so this might not be a very realistic experimental setup, at least in the present configuration. 



minosities required to achieve a useful physics reach point towards ions with small 
charge to minimise space charge and half-lives ~ 1 second to reduce ion losses during 
the acceleration stage whilst maintaining a large number of useful decays per year. 
The most promising candidate ions are ^^Ne and for neutrinos, and ^He and ^Li 
for antineutrinos. A variant on the beta-beam idea is the use of electron capture to 
produce monoenergetic neutrino beams. Electron capture is the process in which an 
atomic electron is captured by a bound proton of the ion A{Z, N) leading to a nuclear 
state of the same atomic number A, but with the exchange of the proton by a neutron 
and the emission of an electron neutrino, 

A{Z,N) + e' ^A{Z-l,N + l) + u,. (7) 

The i dea of using this process in neutrino experiments was independently discussed in 
j^gfg I31J331 j^gf 34|^ jQj^g -^[^j^ Qec- value and long half-life, such as ^^°Sn, were 
proposed to be accelerated to very high boosts with the LHC Baselines of 250 km 
and 600 km were considered with the spectral information coming from the position 
of the events in the detector. Sensitivities comparable to a Neutrino Factory were 
obtained for a single boost. However, in order for electron capture machines to become 
operational, nuclei with shorter half-life are required. The recent discovery of nuclei 
far from the stability line with kinematically accessible super-allowed spin-isospin 
transitions to giant Gamow- Teller resonances (see, for example, Ref. opens up 
such a possibility. The rare-Earth nuclei above ^^^Gd have a short enough half-life 
to allow electr on cap ture processes in the decay ring, in contrast to fully-stripped 
long-lived ions . This was the scenario put forward in Ref. where the use 

of short-lived ions with Qsc-values around 1-4 MeV was proposed. Machines such as 
the SPS, an upgraded SPS and the Tevatron could then be used for the acceleration. 
The ion ^^°Dy, with QEC-value 1.4 MeV, was investigated for the CERN-Frejus (130 
km) and CERN-Canfranc (650 k m) bas elines and different boost factors. It was found 
to have very good physics reach EI132J Owing to the monochromatic nature of the 
beam, multiple boosts are necessary to resolve the intrinsic degeneracy in this case. 

In the following, we demonstrate how the flux for the electron capture/beta-beam 
can be built up by discussing them separately and comparing branching ratios. Let 
the mass difference between the parent and the daughter atoms, AM^ = Ma{Z, N) — 
Ma{Z — 1, + 1), include the mass and the binding energy of an atomic electron 
as well. For electron capture, the maximum kinetic energy release is thus given 
by Qec = AM^^. For /?+-decay, however, the final atom has an excess electron 
since a positron is produced. The maximum kinetic energy release is thus given by 
Q/3+ = AM^ — 2me. Clearly for (AM^^ =) Qec < 2me, electron capture is the 
only allowed process for a proton- rich nucleus. For Qec > 2me, electron capture 
and positron emission compete, their branching ratios dependent on Qec- If decay 
through a emission is also allowed, it is important that this has a relatively low 



Q-value so as not to be the dominant channel For a num ber of useful ion decays 
per year A^^ions; the electron capture neutrino flux is given by I 31 | 32 | 
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for each decay channel. Here, L is the baseline, 7 is the Lorentz boost, E^'^ (= Qec) 
is the neutrino energy in the ion rest frame and E^, is the neutrino energy in the lab 
frame. 

The flux for the /3-spectrum is found in the usual way. In the rest frame of the 
ion, the electron neutrino flux is proportional to 
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Here, E^ 



Qf3+ + rUe = Qec — T^e) IS the total end-point energy of the decay. The 
neutrino flux per solid angle at the detector located at distance L from the source 
after boost 7 is 
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Similarly to the case of electron capture, a neutrino with energy E^-f in the rest frame 
will have a corresponding energy E^, = 2'yEj.f in the laboratory frame along the ^ = 0° 
axis. 

All the known nuclear structure information on the A = 148 and A = 156 nuclides 
has been reviewed in Ref. and Ref. respectively, where the information ob- 
tained in various reaction and decay experiments is presented, together with adopted 
level schemes. Currently, a systematic study of electron capture decays in the region 
of ^"^^Gd, relevant for monoenergetic neutrino beams, is being carried out^^^. Here, 
we consider two nuclides, 7o^Yb and gg^^Tb, that decay through electron capture 
and /3^-decay with similar branching ratios whose lifetimes are not too long or too 
short. Their decays are summarised in Tables [T] and [2J Ytterbium is a nuclide Yb 
with spin-parity O"*", which decays 90% via electron capture plus P^-decay^^ , with 
38% via electron capture and 52% via /3^-decay^^. The remaining 10% goes into 



''The a decay branching ratio is strongly dependent on the QEC-value. For low QeC; the a decay 
probability is sufficiently long as to allow the weak decay modes to be the main channels. 



Decay 


Daughter 


Neutrino energy (MeV) 


BR 


EC 

a 


156 T^™* 

69 -"-'1^ 

156 T^™* 
69 -"-'1^ 

152 p 
68 


2.44 (endpoint) 
3.46 


52% 
38% 
10% 



Table 1: Decay summary for 70^ Yb. The QEC-value for the transition between ground states is 
3.58 MeV and taking into account the excitation energy of the final nuclear state (0.12 MeV), the 
effective Q||^-value is 3.46 MeV. 



Decay 


Daughter 


Neutrino energy (MeV) 


BR 


EC 


148p J* 


2.05 (endpoint) 
3.07 


32% 
68% 



Table 2: Decay summary for 55 ™Tb. The QEC-value for the transition between ground states is 
5.77 MeV and the effective Q'^Q-vahxe to the excited state is 3.07 MeV. 



a-particles and a different final state. This relatively small branching ratio into a's 
helps the nuclide to have a short enough half-life, 26.1 seconds. It is important to 
note that this electron capture-/5"'"-decay transition has only one possible daughter 
state with spin-parity l"*", i.e., it is a Gamow- Teller transition into an excited state of 
Thulium, gg^Tm*. The transition Qsc-value isB QEC-value = 3.58 MeV. However, the 
excitation energy of the final nuclear state (0.12 MeV) needs to be taken into account 
and thus, the effective QEC-value (difference in the total kinetic energies of the system 
after and before the decay) is 3.46 MeV^^. The electron capture energy of ~ 4 MeV 
is well suited to the intermediate-baselines of Europe and the USA with the available 
technology, or those available with future upgrades. On the ot her hand, the gg^^Tb 
isomer with spin-parity Q"*" has a Qec- value of 5.77 MeV^^^^^. Although the decay 
to the ground state of g4'^Gd is highly forbidden, the presence of a Gamow- Teller 
resonance allows the decay into an excited state with effective Q- value 3.07 MeV^^. 
This nuclide is longer lived than yg^Yb (with a half-life of 2.2 minutes) and will require 
slightly higher boosts. It is still well suited to intermediate baselines. However, the 
dominance of the electron capture over the /3"'"-decay channel makes this nuclide less 
desirable. The count rate will be dominated by the single energy of the electron cap- 
ture which provides insufficient information to obtain the g ood sen sitivities aspired to 
by future long baseline experiments. It was shown in Refs. I "^^ ^ "^^ !' that two runs with 
different boosts are necessary for an exclusive or dominant electron capture channel 
to break the intrinsic degeneracy and achieve good CP-violation discovery. Hence, in 
what follows we will study this hybrid approach focusing on ^^^Yb. 



"^QEC'Values are typically calculated between ground states unless stated otherwise. 



Machine 


Tmax 


27maxQ|"c(GeV) 


27maxQ^"+ (GeV) 


SPS 
Upgraded SPS 


166 
369 


1.15 
2.55 


0.81 
1.80 



Table 3: Maximum boosts and neutrino cndpoint energies for ^'^^Yb available for the current SPS 
setup and a proposed 1 GeV upgraded SPS. 



4. Choice of the boost and baseline 

We consider the use of a neutrino beam sourced from boosted ^^^Yb ions directed 
along a single baseline. As described above, both the electron capture and /?"'"-decay 
channels are to an excited state of ^^^Tm with a QEC-value of 3.46 MeV. In order to 
fully exploit the electron capture decay mode, the nuclides cannot be fully stripped; 
at least 16 electrons being left on the ion^^. The maximum boost, 7max; available 
is thus 

-Eacc Z — 16 ,^ 

7max= , , 12 

where mp is the mass of the proton and E^cc is the maximum energy accessible with 
the accelerator. Current and future accelerator facilities would be an ideal production 
environment. In this analysis, we consider the maximum boosts available from the 
current SPS and upgraded SPS (see Table [3]) for the following baselines: 

1. Boost 7 = 166 with current SPS 

• CERN-Frejus (130 km) 

• CERN-Canfranc (650 km) 

2. Boost 7 = 369 with an upgraded SPS 

• CERN-Canfranc (650 km) 

• CERN-Boulby (1050 km) 

With the current magnetic rigidity of the SPS, the electron capture spike can be 
placed on first oscillation for the CERN-Canfranc baseline (650 km) with the beta- 
beam spectrum peaking around the second oscillation maximum (see Fig. [1]). For the 
upgraded SPS, with proton energy at 1 TeV, the first oscillation maximum and most 
of the second oscillation are covered for the CERN-Boulby baseline (1050 km). 





Figure 1: Top panel: appearance probabilities for the CERN-Frejus (130 km) and CERN- 

Canfranc (650 km) baselines. The unoscillated Ve flux in the laboratory frame is shown for ^^^Yb 
given a boost 7 = 166 in arbitrary units. Bottom panel: Ve — > appearance probabilities for the 
CERN-Canfranc (650 km) and CERN-Boulby (1050 km) baselines. The flux from a boost 7 = 369 
is shown in arbitrary units. In both cases, the blue lines correspond to CERN-Canfranc; the red 
being CERN-Frejus (top panel) and CERN-Boubly (bottom panel). The solid lines correspond to 
5 = 0°, dashed 5 = 90° and dotted 5 = -90°. The value sin^ 26*13 = 0.01 was taken for all curves. 



5. Results 

We will compare the physics reach for different experimental Setups, defined by 
the following: 

1. 50 kton detector (LAr or TASD) with 2 x 10^^ ions/yr 

• Setup I: CERN-Frejus (130 km) and 7 = 166 

• Setup II: CERN-Canfranc (650 km) and 7 = 166 

• Setup 111: CERN-Canfranc (650 km) and 7 = 369 

• Setup IV: CERN-Boulby (1050 km) and 7 = 369 

2. 0.5 Mton water-Cerenkov detector with 2 x 10^^ ions/yr 

• Setup III-WC: CERN-Canfranc (650 km) and 7 = 369 

• Setup IV- WC: CERN-Boulby (1050 km) and 7 = 369 

1 and 11 correspond to present SPS energies for the boost, whereas 111 and IV need an 
upgraded SPS with proton energy 1 TeV. The difference between lll(lV) and lll(lV)- 
WC is in the detector. In the first case, we consider a 50 kton detector, like LiAr or 
TASD, with energy reconstruction and neutrino spectral information. In the second 
case, a 0.5 Mton Water Cerenkov detector with neutrino energy from QE events only 
plus inelastic events in a single bin, with 70% efficiency. Whereas in the ffist case 
we use then the total cross section at each energy, for the WC detector we use the 
total cross section at the EC spike, the QE cross section at each energy bin and the 
inelastic cross section for a single bin in the f3~^ spectrum. 

The separation between the energy of the EC spike and the end point energy of 
the beta-spectrum is possible: if Ei^^QE) > 2'yEQ, since i?*™^ > E^{QE), the event 
has to be attributed to the EC flux and, hence, it is not necessary to reconstruct the 
true neutrino energy. 

5.1. Comparing baselines 

For this comparison, we take the combined beta-beam and EC fluxes with ^13 = 1° 
and 6 = 90° for Setups 1 and 11, with different baselines. The result is shown in Fig. [21 
where we see that the longer baseline (650 km) is necessary in order to separate the 
mixing parameter and the CP phase. Violation of CP can be established, only in the 
last case, for a limited range of values of the parameters. 

In fact, the beta-beam channel contributes very little to the overall sensitivity 
of the setup. This is due to the energy dependence of the flux. The smaller flux. 
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Figure 2: 90%, 95% and 99% CL contours for setup I (left panel) and setup II (right panel). The 
parameters 6*13 = 1° and S — 90° have been taken assuming normal mass ordering and 6*23 = 45°. 

combined with the smaller cross section at the energies centred on second oscillation 
maximum, supplies a scarce count rate. The bulk of the sensitivity is due, in this 
case, to the electron capture channel, placed on first oscillation maximum, as seen in 

Fig.m 

5.2. Comparing boosts at the same baseline 

We are going to compare the results for the setups 11 and 111. In going from 7 = 166 
to 7 = 369, the electron capture beam is placed in the tail of first oscillation maximum. 
This gives the beta-beam coverage of the second oscillation maximum and substantial 
portions of the first oscillation maximum. Then the roles of electron capture and beta- 
beam are reversed. The beta-beam now contributes much more to the sensitivity as 
it provides substantial information from the first oscillation maximum and a much 
higher count rate from the second oscillation maximum. 

In Fig. [3] we show the 90%, 95% and 99% CL contours for setups 11 and 111 with 
mixing ^13 = 1" and S = 90°, using the combination of beta-beam and EC fluxes. We 
conclude that the sensitivity is better with the upgraded SPS energy. 

5.3. The virtues of combining energies from beta-beam and EC 

The full power of the combination between the beta-beam spectrum and the EC 
channel is best illustrated in Fig. HI For Setup 111, we give the results with ^^13 = 3° 
and 6 = 90° for the contribution of the beta-beam, for that of the EC channel and 
for the combination. Each of the two techniques separately suffer from a continuum 
of solutions. In fact, the shapes of the allowed regions can be understood by looking 
at the form of the oscillation probability. The power of the combination of the two 
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Figure 3: 90%, 95% and 99% CL contours for setup II (left panel) and setup III (right panel). The 
parameters ^13 = 1° and S — 90° have been taken assuming normal mass ordering and 623 = 45°. 
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Figure 4: 90%, 95% and 99% CL contours for the setup III, simulated for 6I13 = 3° and 6 = 90° 
assuming normal mass ordering and 623 = 45°. The left panel is the contribution of the beta-beam, 
the middle panel is for the electron capture channel with the right panel being the combination. 



channels under the same conditions is in the difference in phase and in amphtude 
between the two fake sinusoidal solutions, selecting a narrow allowed region in the 
parameter space, much more constrained than the two separate techniques. 

The marked difference between the beta-beam alone and the combination with 
the electron capture in this case demonstrate the importance of data from the high 
energies. 

5.4- Disentangling 61^ and 6 

In Fig. [5], we show the physics reach for setup III-WC, with the effects of the 
hierarchy clone solution taken into account. From a comparison of Figs. H] and [5l 
the increase in event rates improves the results substantially, although not as much 
as the size factor between the two detectors. However, owing to the relatively short 
distance, L = 650 km, the mass ordering can be determined only for large values of 
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Figure 5: 90%, 95% and 99% CL contours for setup III-WC with solutions from discrete degeneracies 
included for 6*13 = 1° (left panel) and ^13 = 3° (right panel) for different values of the CP-phase, 
S = -90°, 0°, 90°, 180°. 



the mixing angle ^13. The hierarchy degeneracy worsens the abihty to measure ^13 
and 6 with good precision, especially for negative true values of 6. We conclude that 
a baseline L = 650 km, at least, is needed to disentangle the CP phase. 

6. CP-violation discovery potential 

We now consider the comparison of setups III-WC and IV- WC for the separation 
of the CP phase 6. Boulby provides a longer baseline L = 1050 km than Canfranc 
L = 650 km. This has two contrasting effects on the sensitivity to measure CP- 
violation: i) Sufficient matter effects to resolve the hierarchy degeneracy for small 
values of ^13; ii) It decreases the available statistics. 

From Fig. [6], we see that the smaller count rate results in a poorer resolution for 
the longer baseline. However, the longer baseline allows for a good determination of 
the mass ordering, thus eliminating more degenerate solutions. 

In Fig. [7] we give the CP- violation discovery potential for these two setups. Com- 
paring the two locations of the detector, we notice that the shorter baseline (CERN- 
Canfranc) has a slightly (significantly) better reach for CP-violation at positive (neg- 
ative) values of S than the longer baseline (CERN-Boulby). The longer option, how- 
ever, performs slightly better at negative S if the hierarchy is known to be normal 
and significantly better if the ordering is not determined. This is because the longer 
baseline can identify the neutrino mass hierarchy for these values of 613, therefore 
resolving this degeneracy. 
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Figure 6: 90%, 95% and 99% CL contours for setup III-WC (left panel) and IV-WC (right panel), 
with solutions from discrete degeneracies included for 613 = 1°, for different values of the CP-phase, 
S = -90°, 0°, 90°, 180°. 




Figure 7: CP-violation discovery potential at 99% CL for setup III-WC (left panel) and IV-WC 
(right panel). In each case, we present the results for the beta-beam only (blue dotted lines) and the 
combination with the electron capture result (red solid lines), both without (thin lines) and with 
(thick lines) taking the hierarchy degeneracy into account. 




Sin^(26,3) Sin^(29,3) 




Figure 8: Fraction of S for which the neutrino mass hierarchy can be determined at 99% CL for 
setup III-WC (left upper panel) and IV- WC (right upper panel), with present priors in the known 
parameters. The same fraction of S is presented in the lower panel for setup III-WC with negligible 
errors in the known parameters. In each case, we present the results for the beta-beam only (blue 
dotted lines) and the combination with the electron capture result (red solid lines). 



7. Mass hierarchy determination 

In Fig. [8], we present the results for the neutrino mass hierarchy determination 
for the setups with a 0.5 Mton WC detector (setups III-WC and IV- WC). We do not 
consider the CERN-Frejus cases; the shorter baseline being unable to distinguish the 
type of hierarchy. 

In both cases, the contribution from the beta-beam channel is shown in blue 
dashed lines and the result for the combination with the electron capture channel is 
shown by the red solid lines. As matter effects are more important at high energies, 
we see that the inclusion of the electron capture flux improves the results. 



The CERN-Boulby baseline, with its larger matter effect, represents a much more 
promising setup for the determination of the mass hierarchy. Its resolution would be 
possible for all values of S for sin^ 2^13 ~ few x 10~^, even for the present priors in the 
known parameters. The power of having negligible errors in the known parameters is 
presented for setup III-WC associated with the CERN-Canfranc baseline, showing a 
much higher fraction of 5 for which the neutrino mass hierarchy can be determined. 



8. Conclusions 

We have exploited in this presentation the power of using different neutrino en- 
ergies, in a unique experimental setup, in order to disentangle the CP-odd term in 
the interference between the atmospheric and solar amplitudes for the suppressed 
neutrino oscillation P(z/e ~^ ^n)- This strategy has been implemented by the use of 
a parent ion with two comparable channels of decay: electron capture (EC) and (3'^- 
decay (BB) . The CP phase sensitivity is thus obtained by only using neutrinos, thanks 
to the Energy Dependence of the oscillation probabihty provided by the combination 
of the two EC and BB channels. 

Wc found that the two separate channels EC and BB have a limited overlap of 
the allowed regions in the (^^13, 5) plane, resulting in a good resolution of the intrinsic 
degeneracy. With the SPS upgrade to higher energy [Ep = 1000 GeV), one gets a 
better sensitivity to CP- violation discovery and measurement, iff accompanied by a 
longer baseline of 600 km at least. This aim is the main focus for the third generation 
neutrino oscillation experiments. 

The best E/L for higher sensitivity to the mixing |f/(e3)|, typically the first 
oscillation maximum, is not the same than that for the CP phase determination. 
Like the phase-shift in a interference pattern, the effect of 5 is easier to observe by 
going to the energy region of the second oscillation. This covering is obtained with 
the combination of the two EC and BB channels for the same ion and the same boost. 

The setups 111 and 111-WC, with a baseline L = 650 km, have larger counting 
rates and a better tuning of the beam to the oscillatory patern, resulting in a very 
good ability for disentangling the two parameters (6'i3,5). These setups provide the 
best sensitivity for CP- violation with positive values of S. 

For negative 5's, the type of hierarchy cannot be still resolved for L — 650 km. In 
going to setups IV and IV- WC, with a baseline L = 1050 km, the determination of 
the mass hierarchy is better and one obtains a good reach to CP- violation for negative 
S. 

The general conclusion is that the combination of the two EC and BB channels 
from a single decaying ion and a fixed 7-boost achieves remarkable results for the 
discovery of CP- violation and the measurement of the CP-phase. This is a virtue of 



the different energy dependence of the CP-odd term and the CP-even terms in the 
oscillation probability, and there is no need of performing separate experiments, with 
different systematics and counting rates, for neutrinos and antineutrinos. 
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